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' Abstract. A systematic study of the linear thermal instability of a self-gravitating 

, magnetic molecular cloud is carried out for the case when the unperturbed back- 

ground is subject to local expansion or contraction. We consider the ambipolar 
' diffusion, or ion-neutral friction on the perturbed states. In this way, we obtain 

a non-dimensional characteristic equation that reduces to the prior characteristic 
equation in the non-gravitating stationary background. By parametric manipulation 
of this characteristic equation, we conclude that there are, not only oblate conden- 
sation forming solutions, but also prolate solutions according to local expansion or 
contraction of the background. We obtain the conditions for existence of the Field 
\ lengths that thermal instability in the molecular clouds can occur. If these conditions 

^ . establish, small-scale condensations in the form of spherical, oblate, or prolate may 

' be produced via thermal instability. 

■ 
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1. Introduction 



• With tlie increase in the observational resolution, smaller sizes of den- 

'to '. sity fluctuations have been detected. Direct imaging of CO in nearby 

^ ' clouds revealed substructures with scales of ~ O.Olpc (~ O.OIMq) (Peng 

et al. 1998, Sakamoto & Sunada 2003). Studies of the time variability 
^ \ of absorption lines indicates the presence of small-scale clumps in the 

;h ■ dense gas on scales so small as ~ 5 x lO^^pc (~ l^AU), with masses of 

■ ~ 5 X 10-^Mq (Moore & Marscher 1995, Rollinde et al. 2003). 

As a general rule, neither subparsec nor AU-scale condensations 
in molecular clouds are spherical (Ryden 1996). Jones & Basu (2002) 
have recently decipher intrinsic three-dimensional shape distributions 
of molecular clouds, cloud cores, Bok globules, and condensations. They 
find out that molecular clouds mapped in ^^CO are intrinsically triaxial 
but more nearly prolate than oblate, while the smaller cloud cores, Bok 
globules, and small-scale condensations are also intrinsically triaxial 
but more nearly oblate than prolate. 
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Small-scale condensations appear to be immediate precursor of large- 
scale clumps (dense cores with significant jeans mass) via merging 
and collisions; they constitute the initial conditions for star forma- 
tion. Therefore, the understanding of the origin and merging of these 
small-scale condensations is of fundamental importance for a consistent 
theory of star formation and galactic evolution. 

The origin and shape of these small-scale condensations is a dis- 
putable issue. Anisotropic heating and fragmentation via gravitational 
collapse is an important reason for oblate/prolate large-scale clumps 
with significant Jeans mass (e.g. Nelson &: Langcr 1997, Indebetouw 
&; Zweibel 2000, Hartmann 2002). The above scenario is not correct 
for small-scale condensations, because they have low gas density and 
small sizes, thus, their masses are significantly smaller than their cor- 
responding Jeans mass. According to this feature, the only remaining 
responsible parameters may be turbulence and/or thermal instability. 

Gammie et al. (2003) have recently studied the effect of turbulence 
in three dimensional analogs of clumps using a set of self consistent, 
time-dependent, numerical models of molecular clouds. The models 
follow the decay of initially supersonic turbulence in an isothermal, self- 
gravitating, magnetized fluid. They have concluded that nearly 90% of 
the clumps are formed in prolate and 10% of them are oblate. 

In molecular clouds, the dispersion velocity inferred from molecu- 
lar line width is often larger than the gas sound speed inferred from 
transition temperatures (Solomon et al 1987). Magnetohydro dynamic 
turbulence may be responsible for the stirring of these clouds (Arons 
Sz Max 1975). Because of these turbulent motions, molecular clouds 
must be transient structures, and arc probably dispersed after not much 
more than ~ lO^yr (Larson 1981). Since cooling time-scale of molecular 
clouds is approximately ~ 10^ — lO^yr (Gilden 1984), thermal instability 
may be a coordinated trigger mechanism to form condensations. Turbu- 
lence, in the second stage, can deform these small-scale condensations 
in shape and orient them relative to the background magnetic fields. 

Observations and theoretical studies establish that magnetic fields 
play an important role in shaping the structure and dynamics of molec- 
ular clouds and their substructures (e.g. Basu 2000, Fiege & Pudritz 
2000, Hennebelle 2003). The relative alignment of the projected mag- 
netic field with the projected minor axis of the condensations is an 
important diagnostic. 

In conformity with the above explanation, Nejad-Asghar &: Ghan- 
bari (2003 hereafter NG) interested to investigate the effect of ambipo- 
lar diffusion on the thermal instability and formation of small-scale 
condensations in the magnetic molecular clouds. They concluded that 
there are solutions where the thermal instability allows compression 



nejadasghar.tex; 5/02/2008; 10:23; p. 2 



Small-Scale Condensation Formation 



3 



along the magnetic field but not perpendicular to it. NG inferred that 
this aspect might be evidence in formation of the observed oblate 
small-scale condensations in magnetic molecular clouds. 

In this paper we want to testify and develop the work of NG by 
including self-gravity and local background contraction/expansion. We 
present the basic equations, background evolution, and the linearized 
equations in section 2. Section 3 deals with exponential growth rate and 
parametric solutions that culminates in formation of oblate, prolate, 
and spherical condensations. Section 4 allocates to a conclusion and 
some future prospects. 



2. The Equations of The Problem 

The basic equations, including self-gravity and ambipolar diffusion, 
are given first in general (§ 2.1) and then specialized for the homo- 
geneous contracting/expanding molecular cloud (§ 2.2) and for small 
perturbations to that medium (§ 2.3). 

2.1. Equations 

In principle, the ion velocity vi and the neutral velocity Vn in molecular 
clouds, should be determined by solving separate fluid equations for 
these species (Draine 1986), including their coupling by collision pro- 
cesses. But, in the time-scale of cooling considered here, (10^ — lO^yr, 
Gilden 1984), two fluids of ion and neutral are well coupled together, 
and we can use the basic equations as follows (Shu 1992) 

| + pV.. = (1) 

+ Vp + V(|J) - (B • V)^ + = (2) 
' ^ ^'^P + p^-V.{KVT) = ^ (3) 



^ — 1 dt ^ — 1 p dt 
(7B B 

— + B(V-t;)-(B.V)t; = Vx{^-— j^x[Bx(VxB)]} (4) 

= 47rGp (5) 
p - -pT = (6) 

where variables and parameters have their usual meanings and r] k, 
2.46 X 10^^cm^.5~^.s~^ is the collision drag in molecular clouds (see 
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Figure 1. Logarithm of the coohng rate, A(no, To) = AonoT^, versus number density 
in molecular clouds, no{cm~^). The values of 5 and (3 are shown in this figure. 

McDaniel & Mason 1973). We use the relation pi = ep'^ (e ^ 1.83 x 
10~^''cm~^/^.5^/^, V = 1/2) between ion and neutral densities (Ume- 
bayasi k. Nakano 1980), and in a good approximation we choose p = 

Pn + Pi^ Pn- 

Q{p,T) = A(p, T) — Ftot is the net cooling function [erg.s^^ .g^^), 
where Ftot is the total heating rate and A{p, T) is the cooling rate which 
can be written as (Goldsmith & Langer 1978, Neufeld et al. 1995) 

A(p,r)=Aop''r'^ (7) 

where Aq. S, and arc constants. The range of (3 is 1.4 to 2.9. The 
constant 5 is greater than zero for optically thin case and less than 
zero for optically thick case (see Fig. 1). 

Models of the molecular clouds identify several different heating 
mechanisms. In this paper, we consider the heating rates of cosmic rays, 
H2 formation, H2 dissociation, grain photoelectrons, and collisions with 
warm dust, as a constant Tq (Glassgold & Langer 1974, Goldsmith &; 
Langer 1978). The heating of the gas by magnetic ion-neutral slip is 
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discussed in detail by Scalo(1977); a simple estimate of this heating 
rate is 

Tad = vep''^ (8) 
where va is the drift velocity of ions 

In order of magnitude, if kBq changes on a typical scale of A, then 

and ambipolar diffusion heating rate is given by 

Tad = T'op-^^^'^) (11) 

where Fq is defined as 

^°-167rVA2- ^ ^ 

The gravitational heating rate is found by setting the rate of compres- 
sional/expansional work per particle, pd{n~^)/dt, equal to the rate of 
change of gravitational energy per particle, [d{PEtot)/dt]/{nV), where 
PEtot is the gravitational potential energy of the volume V. For a 
uniform sphere of radius A, we find 

rgrav = r'y/^ (13) 

where Fg is defined as 

(fvrG)^ 2 .... 

where d(r) is the contraction/expansion parameter rate (see §2.2). 
2.2. Background Evolution 

As a basis for the small-perturbation analysis, we assume a homoge- 
neous background which is expanding/contracting uniformly, so that 
the unperturbed quantities only depend on time. The background quan- 
tities will be denoted with the subscript 0. The expansion is given 

by 

r = a{t)x (15) 
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where r is the Eulerian coordinate, x is the Lagrangian coordinate and 
a{t) is the expansion/contraction parameter. Using equation (15), the 
unperturbed velocity field is given by 

VQ{r,t) = ^^^r. (16) 
a 

for the background evolution, the basic equations (l)-(6) reduce to 

Po(i) = Po(t = 0)a(t)-3, pn{t) = po{t = 0)a(t)-3T 
ro(t) = ro(i = 0)a(t)-3(T-i), Bo{t) = Bo{t = 0)a{t)-^ (17) 

where a{t) follows the differential equation 

airfdir) = -1 (18) 

where the dot over the symbol indicates the derivative respect to a 
non-dimensional variable r = [^irGpo{t = 0)]^/^t. 

2.3. Linearized Equations 

Density fluctuation ratios in the molecular substructures is in the order 
of ~ 10 (Falgarone et al. 1992, Pan et al. 2001). Therefore, the linear 
regime of the thermal instability might lead to some significant results 
for small-scale condensation formation. 

To obtain a linearized system of equations, we split each variable into 
unperturbed and perturbed components, indicating the latter with a 
subscript 1. Eulerian del operator is applied to the equations and all 
equations arc rewritten in terms of the Lagrangian coordinate x. The 
resulting linear system has coefficients which depend on t but not on x . 
We then carry out a spatial Fourier analysis, with Fourier components 
proportional to exp{ik • a;), so that k is the Lagrangian wave vector. 

Simplifying the resulted linear system by repeated use of the back- 
ground equations (17); we obtain 

^(Pl)+ik.vi=0 (19) 
at Po 



dvi da/at . Cs k pi .Bq-Bi-- 
at a ^Ts k Po 47r/9o 



-.^B,-z^(^)=0 (20) 
47rpo k r^ po 



9 



(tl) + . + _ + + ^ = 21 

dt Po Tcp TK Po TcT Tcp Tr PO 
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+ tBo{k ■ vi) - i{k ■ Bo)t;i + ^Bi 



dt 



+fe X { 



Bn 



47rr?epo 



i+i 



- X [Bo X (fc X Bi)]}= 



(22) 



where Cs = vTPo/po and k = k/a are, respectively, the adiabatic back- 
ground sound speed and the Eulerian wave vector. The other symbols 
have the following definitions: 



T. = 



kc.i 



Tn = 



TK = 



Rpo 



2 ' 



/x(7- l)po{dn/dp)T' 



/x(7 - l)Kk: 
R 



(23) 



that are the characteristic time-scale of sound waves, self-gravity per- 
turbation waves, thermal conduction, isothermal and isobaric differen- 
tial cooling, respectively. 

We use the coordinate system Ux, Uy, and Uz as specified by NG. 
Equations (20) and (22) may be used to uncouple viy- the perturbed 
velocity in the plane perpendicular to both Bq and k- from the rest 
of the problem. With the choice of exponential perturbation {e^^), dis- 
turbances perpendicular to the (Bq — fe)-plane, have a solution which 
displays existence or non-existence of the Alfven waves. Amplitude of 
the Alfven waves are damped via expansion of the medium and/or 
with ion-neutral friction, while, it must grow with injection of energy 
in contracting medium. 

The motion in the other modes are constrained to the xz-plane, and 
are governed by the matrix equation. 



y(i) = AY 



(24) 



where y is a 5 x 1 matrix as 



Y = 



Pi/po 
Pi/Po 
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and y(^) is its first time derivative. The 5x5 matrix of the coefficients, 
A, is defined as 



A 













1 


1 


1 




TcT 




TK 


ic.,T.,asm 




ic.,asin 


e 


'g 








iCsTsacos 


9 


icsocos 


e 


'g 



















CsTga 

i-y sin 6 

CgTsa 








V 

where 6 is the angle between k and Bq, and 
1 



\/47rpo 



TAD = 



_ I cost) 

CsTsa 

il cos 6 
CgTsa 











iCsTsa 
1-'^ 



\ 



1 



(25) 



are the characteristic time-scales of the Alfven waves and ambipolar 
diffusion, respectively. 



3. Exponential Growth Rate 

The standard exponential growth rate provides the following formal 
solution for all the perturbations: 

yi(t)=yiit = 0)exp{ht), (26) 

where real{h) represents the growth/decay rate. Inserting the back- 
ground evolution (eq. [17]) and exponential growth form (eq. [26]), into 
equation (24); we obtain 

= hY + CY (27) 
where C is a diagonal matrix as 

C = —dm5[3, 37, 1,1,2] (28) 

where | Tg |= |^^°^^| represent contraction/expansion time-scale. Exis- 
tence of solution for equation (27), needs the following condition: 

Det[hl + C-A]={) (29) 

where / is the unitary matrix. By introducing the non-dimensional 
quantities 

y = hTs, (Tp = — , or = — + — , 

TcT Tcp TK 

a ^ {^f, D ^ = (^)^ Ec = ^, (30) 

^TAL K^TAD T„ Te 
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we find a five-degree finear characteristic equation that without self- 
gravity and expansion/contraction of the background (Gg = 0, Ef> = 0), 
reduces to the equation (22) of NG. We use the Laguerre method to 
find the roots of this characteristic equation. 



3.1. The Sound Domain 



When the sound period, Tg, is much smaller than the other character- 
istic time-scales, we can neglect the effect of a, D, Gg, and E^. The 
characteristic equation reduces to a three-dcgrcc Uncar equation. There 
are three solutions: two sound waves and one condensation mode. The 
stable region for this case is shown in the Fig. 1 of NG. 

For the temperatures and densities that thermal instability is desta- 
bilizing, there exist a critical length (Field length) defined as the maxi- 
mum wavelength which thermal conduction can suppress the instability. 
Inserting the defined time-scales (eq. [23]), into the definitions of ap and 
or; we obtain 

= ^i3,7^^ A(no, n)i6 + 2.5^ - 1.5x) (31) 



^T = ^^^^Hno,Tomi + {^f] (32) 



where ^ and % are 



Fad 

^ = A(no,ro)' ^= A(no,ro) ^ ^ 



and Ao is a defined wavelength as follows: 



The values of the Ao(pc) for typical data in the molecular clouds is 
given in Fig. 2. According to this figure, we choose wavelengths in the 
range of A(pc) ~ 10~^ — 10~^pc, that are interesting in formation of the 
small-scale condensations. 

We can now define a generalized Field length for two cases as follows: 

— 6 + 2.5^ — 1.5% > /?, which ap > 0, that is upwards of the ax — CTp 
plane: 

A» = ^ ^° (35) 

^ / ■5+2.5|-1.5x ~ ^ ^ 
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2.0 2.5 3.0 3.5 4.0 4.5 5.0 



log(no) 

Figure 2. Logarithm of the Field length (at unit of parsec) in the typical molecular 
clouds versus number density, no{cm~^). 



S + 2.5C-1.5X < 
o'T — CTp plane: 



-^0 



which fjp < 0, that is downwards of the 



3 i.5x-<5- 
2 



2.5? 



(36) 



If (5 + 2.5^ — 1.5x sets between — |/? and /?, the Field length can not 
be defined, thus, the medium is stable for all wavelengths. Otherwise, 
for wavelengths greater than the defined Field lengths (eq. [35] and 
[36]), the medium is unstable. This case is shown in Fig. 3 for two 
temperatures of lOi^ and WOK. 

Before proceeding any further, we must have a real physical feeling 
of defined non-dimensional parameters: a, D, E^, and Gg. We consider 
some typical magnetic molecular clouds with density between lO^cm"'^ 
to lO^cm"^, temperatures in the range of Tq ~ 10— 100-R', and magnetic 
field strength Bq 10/iG (Myers & Goodman 1988, Crutcher 1999). 
The typical values of a, D, Ee, and Gg are shown in Fig. 4. 



3.2. The Magnetic Domain 

Wherever the Alfven period and the ambipolar diffusion time-scale 
are important in the magnetic molecular clouds, we can not ignore 
the effect of a and D in the characteristic equation. This case has 
recently been investigated by NG, without self-gravity and local expan- 
sion/contracion of the background {Gg ~ 0, -Eg ~ 0). They conclude 
that there are solutions where the thermal instability allows compres- 
sion along the magnetic field but not perpendicular to it (see Fig. 2 
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2.0 2.5 3.0 3.5 4.0 4.5 5.0 2.0 2.5 3.0 3.5 4.0 4.5 5.0 



log(r\,) log(r^ 
(a) (b) 

Figure 3. The stable, isobaric instability, and isentropic instability of the medium 
which can occur according to the different values of ^ — 0.6x, for (a) To = IQK and 
(b) To = imK. 
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Figure 5. Regions of stability(x), spherical instability(»), and oblate instability (—) 
for the case of Gg = 0.1 and Ee = 0, for a typical value of a = 5.0 and D = 1.0. 
The decreased stability region in the isobaric instability criterion (line OA), occurs 
by the increasing of internal pressure via existence of self-gravity. 



and Fig. 4 of NG). Maximum effect of the magnetic domain for the 
oblate condensation formation is occurred when ambipolar diffusion 
time-scale is nearly equal to the period of sound waves D « 1.0 (i.e. 

3.3. The Self-Gravity and Expansion/Contraction Domain 

In this subsection we investigate the effect of the self-gravity and the 
expansion/contraction of background. If we consider the effect of the 
self gravity without expansion/contraction of the background, the iso- 
baric instability criterion (line OA of Fig. 1 of NG) is modified via 
bringing this line downwards. This case is shown in Fig. 5, for a typical 
value of a = 5.0 and D = 1.0. Physically, this means that self-gravity 
causes to increase the internal pressure, thus, isobaric instability must 
occur at a decreased ap for each or- 
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Figure 6. Regions of stability ( x), spherical instability(»), and oblate instability(— ) 
in the expanding background, for a typical value of a = 5.0 and D = 1.0, with (a) 
Be = 0.01 and Gg = 0.1, and (b) = 0.1 and Gg = 0. 



If the background is expanding {E^, > 0), its expansion energy causes 
to stabilize the medium. This case is shown in Fig. 6 for a typical value 
of CK = 5.0 and D = 1.0. In the isentropic instability criterion (line 
OB), expansion energy causes to stabilize the medium in the direction 
of the magnetic field and perpendicular to it. On the other hand, in 
the isobaric instability criterion (line OA), it only causes to stabilize 
the medium in perpendicular to the magnetic field, corresponding to 
decreased pressure via ion-neutral friction. 

For contracting background (E^ < 0), contraction energy injected to 
the medium, thus, its stability is decreased and converted to a prolate 
instability. Diffusion of neutrals relative to the freezed ions in the per- 
pendicular direction of the magnetic field is the reason of this prolate 
instability. This case is shown in Fig. 7 for a typical value of a = 5.0 
and D = 1.0. When the parameters of a magnetic molecular cloud set, 
locally, in this region of cjt — o"p plane, prolate condensation may be 
produced via thermal instability. 
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4. Summary And Prospects 



In this paper we perform linear analysis of thermal instability in a lo- 
cally uniform expanding/contracting magnetic molecular clouds which, 
in the perturbed state, is undergoing ambipolar diffusion. Thermal 
conduction and self-gravity have also been included as fundamental 
ingredients. The small-perturbation problem yields a system of ordi- 
nary differential equations with five independent solutions. Wc choose 
an exponential growth rate which convert the system of ordinary dif- 
ferential equations into a five-degree complete characteristic equation. 
If we neglect the self-gravity and expansion/contraction of the back- 
ground, the characteristic equation reduces to the prior results of NG. 
We have used the Laguerre method to find the roots of this complete 
characteristic equation. 

In sound domain, two of the solutions have the character of oscil- 
latory modes (sound waves) and the third one is a non-oscillatory (or 
condensation) solution. We adopt a parametric net cooling function and 
find for perturbations with wavelengths greater than the Field length, 
thermal instability causes the medium to condense. Fig. 3 shows the 
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condition of instability in the molecular clouds which their cooling rates 
are presented in Fig. 1. 

We choose a wide range of density and temperature in the molecular 
clouds with typical magnetic field strength Bq « lO/xG. Interesting 
wavelengths in the problem are around the Field length which is shown 
in Fig. 2. According to this figure, we consider wavelengths around 
10"'* — 10~^pc for small-scale condensations. The typical values of a, 
D, Eg, and Gg are shown in Fig. 4. 

In magnetic domain, without self-gravity and expansion/contraction 
of the background, there arc solutions where the thermal instability 
allows compression along the magnetic field but not perpendicular to 
it. Maximum cases of these oblate condensation solutions are occurred 
when the ambipolar diffusion time-scale equals to the period of sound 
waves {D ~ 1). 

Fig. 5 to Fig. 7 show results that come into existence by consid- 
ering of self-gravity and local expansion/contraction of the magnetic 
molecular cloud. We deduce that the self-gravity causes to increase the 
internal pressure, thus, isobaric instability must occurs at a decreased 
ap for each ctt- Therefore, instability of the medium is increased. In 
the expanding background, expansion energy in the isentropic insta- 
bility, causes to stabilize the medium in the direction of the magnetic 
field and perpendicular to it, while, in the isobaric instability it only 
stabilize the medium in the perpendicular direction. In the contracting 
background, stability of the medium is decreased and converted to the 
prolate instability via injection of contraction energy and diffusion of 
free neutrals relative to freezed ions. 

In this paper we conclude that linear thermal instability can produce 
small-scale condensations in spherical, oblate, or prolate. We try to 
analyze a rather involved problem linearly, because, before this process 
overcome, turbulence causes to interaction and merging of these newly 
formed condensations. Physically, we expect that merging of these 
small-scale condensations culminate in the large-scale clumps that are 
star bearing regions in our world. Authors are now preparing a com- 
plete simulated turbulent magnetic molecular cloud with condensations 
produced by thermal instability. It would be interested to investigate 
the effect of interaction, merger, and coagulation of these small-scale 
condensations with smoothed particle hydrodynamics (SPH) method. 
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